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VISIBLE  SURFACE  ALGORITHMS 


FOR  QUADRIC  PATCHES 


ABSTRACT 


This  paper  describes  two  algorithms  which  find  the  visible  por¬ 
tions  of  surfaces  in  a  picture  of  a  cluster  of  three-dimensional  quadric 
patches.  A  quadric  patch  is  a  portion  of  quadric  surface  defined  by 
a  quadrat" c  equation  and  ^y  zero,  one  or  several  quadratic  inequalities. 
The  picture  is  cut  by  parallel  planes  called  scan  planes;  the  visibility 
problem  is  solved  in  one  scan  plane  at  a  time  by  making  "a  good  guess" 
as  to  what  is  visible  according  to  the  visible  portions  found  in  the 
previous  scan  plane. 

The  algorithm  for  intersecting  patches  works  in  a  time  roughly 
proportional  to  the  number  of  patches  involved  (and  not  to  the  square 
of  this  number  as  with  some  previous  algorithms) . 
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INTRODUCTION 


N™  StUdieS  [2'5)  —  ^oted  to  the  problem  of  finding 

visible  portions  of  surfaces  when  the  surfaces  were  defined  as  being 

planar  Polygons.  Watkins  ,2)  desoribes  one  of  these  algorithm  „hich 

“  quite  fast  and  requires  little  ament  of  working  storage. 

Pew  successful  studies  have  been  made,  however,  for  non-planar 
surfaces.  BE  VISION  [1)  finds  the  hidde!n  Unes  „hm  ^  ^ 

defined  by  quadratic  equations  and  by  quadratic  inequalities.  Another 
approach  was  taken  at  the  University  of  Utah  by  Henri  Oouraud  [3J  who 
ppioximates  Coons'  rational  cubic  surface  patches  [4]  with  twisted 

polygons  and  then  uses  Watkins'  alanvithm  • 

algorithm  in  order  to  find  the  visible 

of  surfaces  Tn  j  i 

aces,  in  Gouraud's  method,  the  shading  of  a  point  of 
a  surface  is  obtained  by  interpolating  linearly  the  exact  shading  value 
at  the  four  corners  of  the  surrounding  polygon,  thus  obtaining  ,  contin_ 
uous  shading  function  all  over  a  given  surface. 

The  author  appreciates  the  merits  of  these  previous  approaches,  but 
was  looking  for  a  method  yielding  the  exact  visible  contour,  intersec¬ 
tion  lines,  and  shading  functions.  To  find  this,  he  was  forced  to  limit 
bis  investigations  to  the  same  type  of  quadratic  surfaces  which  were 

d  by  Weiss  [1J  m  order  tc  limit  the  computational  complexity  of  the 

problem  to  the  solution  of  quartic  (4th  \ 

quartic  (4th  degree)  polynomial  equations. 

The  methods  described  in  the  Daae*  fv_n_  • 

Pages  following  are  related  to  the 

Watkins  «  algorithm  (2) .  They  offer  two  kind*  r  f  - 

-  kinds  of  advantages  over  the 

g  im  of  bl  VISION  [1J.  First,  they  yield  visible  surfaces  instead 
of  systematically  eliminating  hidden  lines.  Thus  they  permit  the 
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generation  of  shaded  pictures  with  only  the  further  effort  of  computing 
the  shading  value  along  the  visible  portions  of  the  surfaces.  Secondly, 
they  work  in  an  amount  of  time  which  is  roughly  proportional  to  the 
number  of  surfaces  fed  into  the  algorithm,  rather  than  the  square  of 


this  number. 
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DEFINITION  OF  QUADRIC  PATCHES 

Quadric  patches  are  defined  in  the  same  way  as  was  used  by  Weiss  [1] 
A  quadric  patch  is  a  set  of  points  obtained  by  setting  a  quadratic  expres 
sion  equal  to  zero  and  by  constraining  zero,  one  or  more  other  quadratic 
expressions  to  be  either  positive  or  negative. 

A  quadratic  expression  is  given  by : 

A(x,y,z)  =  a  x2  +  a  y2  +  a  z2  +  a  vz  +  a  zx  +  a  xy  +  a  x  +  a  y 

2  3  *♦  5  6  7  8 

+  a  z  +  a 
9  10 

The  quadratic  inequalities  are  used,  in  particular,  to  define  the 
boundaries  of  the  quadric  patches.  For  example,  suppose  that  we  want 
to  define  an  infinitely  long  cylinder  in  the  z  direction,  centered  at 
(0,0)  and  of  radius  1;  its  equation  would  be: 

(1)  x2  +  y2  -1  =  o 

If  we  were  interested  in  the  portion  of  the  cylinder  truncated  by 
the  planes  z  »  1  and  z  =  -1,  we  would  impose  the  constraints 

(2)  z  +  1  <_  0 

(3)  z  -  1  >_  0 

If,  moreover,  we  want  to  bore  a  hole  through  this  cylinder,  we 
may  impose  the  additional  constraint: 

(4)  x2  +  z2  -  .25  >_  0 

Equations  (1)  to  (4)  define  the  resulting  patch,  shown  in  Figure  1. 


SCAN  LINE  TECHNIQUES 


It  is  necessary  to  solve  the  visible  surface  problem  only  to  a 
certain  limit  of  resolution  which  is  determined  either  by  the  resolution 
of  the  human  eye ,  or  by  the  precision  of  the  scope  on  which  the  picture 
is  to  be  displayed.  This  leads  to  the  idea  of  cutting  the  picture  by  p 
horizontal  lines,  and  computing  the  shading  function  at  p  points  on  each 
line.  The  time  taken  by  the  visible  surface  algorithm  is  proportional 
to  p,  while  the  time  taken  for  shading  is  proportional  to  p2.  The  time 
taken  for  the  exact  computation  of  the  shading  of  quadrics  is  much 
larger  than  for  planar  polygons:  two  minutes  was  typical  on  PDP-10  for 
p  =  500.  The  value  of  p  is  usually  chosen  between  250  and  1000. 

w/e  assume  that  scalings,  translations,  rotations,  and  a  perspective 
transformation  have  previously  been  applied  to  the  objects  such  that  the 
observer  is  looking  along  the  z  axis  toward  smaller  values  of  z  (Figure 
2) .  Lines  of  sight  are  parallel  to  z.  We  consider  in  particular  p 
equidistant  horizontal  planes  (parallel  to  x,z)  called  scan  planes,  which 
intersect  the  transformed  objects.  These  scan  planes  are  projected  onto 
the  image  plane  as  p  equidistant  scan  lines. 

We  solve  the  visible-surface  problem  for  a  given  scan  line  by  cutting 
the  object  surfaces  by  the  scan  plane.  The  intersections  of  quadric 
patches  with  this  plane  are  portions  of  conics.  The  problem  is  then 

the  current  scan  plane,  and  the  values  of  the  shading  function 
are  computed  along  the  scan  line  for  the  visible  portions  of  conics. 

This  method  of  computing  visibility  is  made  fast  because  we  keep  track 
of  the  visible  portions  in  the  previous  scan  plane,  thus  reducing 
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considerably  the  number  of  comparisons  of  portions  of  conics  necessary 
in  order  to  determine  which  one  is  in  front  of  the  others. 

We  are  now  going  to  examine  two  visible  surface  algorithms:  one 
for  non-intersecting  quadric  patches,  and  another  one  for  possibJy 
intersecting  quadric  patches.  If  patches  are  known  not  to  intersect 
one  another,  one  can  cons it er ably  decrease  the  computation  time  for 
visible  surface  detection.  We  shall  examine  such  an  algorithm  first. 

ever,  the  algorithm  for  intersecting  patches  is  more  interesting  and 
involves  some  cleverness  in  avoiding  unnecessary  computations. 
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The  intersection  of  .  quadric  by  a  sen  plane  is  a  conic,  perhaps 
degenerate,  for  the  observer  situated  at  infinity  in  the  .  direction 
a  conic  has  two  portions,  the  closest  Portion  and  the  furthest  portion 

"  *  “"iC  haS  ““  -«*“■  **  the  same  given  x  coordinate 
value,  the  one  with  the  greatest  z  value  belongs  to  the  closest  portion, 
the  other  one  to  the  furthest  portion. 

Spceial  points  of  a  conic  are  of  two  hinds,  apparent  contour  point 

<X‘  “  ^  3'  Uh"e  — -  *"  Perpendicular  to  the  z  axis 

•ad  points  of  intersection  of  the  conic  with  its  boundaries  ,x  and  x  , 

These  points  are  clipped  outside  the  ,0.1,  interval  .points  ouLide  this 
interval  are  eliminated  as  being  off  the  screen,,  they  are  then  ordered 
according  to  increasing  x  values,  for  each  interval  ^ 

visible  portion  of  the  patch  is  —red.  Pot  instant  in^e  3 
these  intervals  are: 


Interval 


lx  -  0,  x  ] 
0  1 

lx  ,  x  ] 

1  2 

lx  ,  x  ] 

2  3 

lx  ,  x  ] 

3  4 

lx  ,  x  =1] 

4  5 


Visibility 
0  (not  visible) 

1  (closest  portion) 

2  (furthest  portion) 
I  (closest  portion) 

0  (no.  visible) 


Ih  the  first  phase  of  the  algorithm,  an  explicit  representation  of 
the  conic  z  -  f.x,  is  also  obtained,  it  win  be  used  to  do  the  depth  coz 
the  next  section  and  later  to  obtain  the  shading  values. 
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THE  ALGORITHM  FOR  NON- INTERSECTING  PATCHES 

We  have  determined,  in  a  scan  plane,  the  visibility  of  each  f  atch 
for  the  observer  if  all  other  patches  were  removed.  Now  wc  must  synthe¬ 
size  all  this  information. 

The  first  step  is  to  merge  together  the  special  points  of  all  the 

patches  into  one  large  list  which  is  ordered  by  increasing  x  coordinates. 

This  list  will  then  be  scanned.  For  each  interval  [x.,  x.  1  we  shall 

1  l+l 

determine  which  is  the  visible  patch,  if  any.  Of  course  at  that  moment 
the  visibility  in  the  [x^_^,  x^]  interval  has  already  been  determined, 
except  for  i  =  0.  Thuss 

if  i  E  0  or  if  x,  is  a  special  point  for  the  patch  which  was 

visible  in  the  [x^,  x^  interval,  we  have  to  compare  all  patches  in 

the  [x^,  x^+^J  interval.  This  is  achieved  by  computing  the  depth  of  each 

patch  (z  coordinate)  in  the  interior  of  the  interval,  for  instance  at 

the  midpoint  x  *=  i  i+1  .  Those  depths  are  then  sorted,  and  the 

2 

patch  with  the  greatest  z  is  declared  visible. 

if  x^  is  a  specia .  point  for  patch  P  where  a  new  portion  of  patch 
P  might  become  visible,  and  if  pc'.tch  P’  w.,a  found  to  be  visible  in  the 
interval  [x^_^,  x^]  with  P'  /  P,  ther.  P  and  P'  must  be  compared,  but  not 
to  any  other  patch.  Of  P  and  P’ ,  the  one  closest  to  the  observer  is 
declared  visible.  If  no  patch  was  visible  in  the  previous  interval,  then 
patch  P  automatically  becomes  visible. 

This  ai.gorithrr  can  be  slightly  improved  by  observing  if  a  patch  A 
was  in  front  of  another  patch  A'  in  an  interval  [xi  ^  x^  and  if  xi  is 
not  a  special  point  for  either  of  these  patches,  the  patch  A  will  still 
cover  A'  in  interval  [x^  +  x^+^] .  Thus  we  can  avoid  unnecessary  depth 


co^ .tatxon,  by  maintaining  deeth;;^  ^  o£  _  Qf  ^  ^ 

WhUe  ““  "  SC""nin9  (i"“—  «>  •  «—  a  special  point  is  found 

£“  3  P3tCh'  PatCl'  ^  dal<!“d  *•  ^-ordered  lists,  except 

”  1£.1S  0,6  “hiCh  WU1  te  —  *"  «-  —  X  interval,  in  which 

case  it  is  placed  in  front  of  all  h0p  mu 

all  lasts .  The  usefulness  of  these  lists 

arises  when  we  arrive  at  a  point  x  whinh 

i  h  h  18  a  sPecial  point  for  the  patch 

“h£Ch  “aS  ^  0,6  ‘■W  V  -us  requiring  all  patches  to 

compared.  However.  patches  which  are  below  the  top  of  any  _  of  our 

depth-ordered  lists  will  certainly  not  he  seen  because  they  are  at  least 

hidden  by  the  patch  which  is  at  the  top  of  that  list  mx 

UP  ot  that  list.  Those  patches  which 

cannot  be  elindnated  in  this  way  will  need  to  be  compared  and  will  be 
depth-sorted  into  a  next  list.  A  depth-ordered  list  is  deleted  whenever 
U  as  reduced  to  only  one  element.  Because  a  given  patch  may  exist  in 
several  of  these  lists  at  a  time,  it  1.  desirable  to  have  all  occurrences 
of  the  same  patch  themselves  list  United  in  order  to  be  able  to  delete 
-cm  rapidly  if  we  arrive  at  a  special  point  for  this  patch. 


1 

A  VISIBLE  SURFACE  ALGORITHM  FOR  INTERSECTING  P  VTCHfc'S 

If  patches  may  intersect,  we  cannot  compare  them  on  an  interval  by 
computing  their  depths  at  some  arbitrary  point  in  this  interval.  The 
correct  way  to  compare  two  patches  on  an  interve 1  is : 

1.  To  compute  their  intersections. 

2.  If  there  is  any  intersection  in  the  interval,  to 
subdivide  the  interval  (creating  n  +  1  subintervals 
if  there  are  n  intersections)  so  that  the  patches  do 
not  intersect  on  each  interval. 

3.  To  compare  depths  at  any  arbitrary  point  within  each 
interval  in  which  the  patches  are  now  known  not  to 
intersect. 

If  n  patches  are  involved,  this  method  would  require  about  n^  /4 
comparisons  in  order  to  determine  which  patch  is  in  front  of  the  others 
oit  an  interval .  The  computation  of  the  intersections  of  two  conics  is 
an  expensive  operation  which  requires  around  3  milliseconds  of  computation 
time  on  the  PDP-10  (on  which  these  algorithms  were  implemented  by  the 
author).  Most  of  this  time  is  devoted  to  solving  a  quartic  equation.  Thus 
it  becomes  critical  to  reduce  the  number  of  comparisons  of  quadrics  and  of 
computations  of  intersections.  This  car  be  done  in  several  ways: 

1.  When  the  intersection  of  two  conics  has  been  determine  t , 
they  must  be  stored  somewhere  in  order  to  avoid  havin 
to  recompute  the  same  intersections  on  the  same  scan 
line.  If  the  numbers  k  and  k'  identify  two  conics,  the 
intersections  of  k  and  k'  are  stores  in  a  hashcoded  array. 

The  access  function  in  this  array  (or  hash  function) 

H (k ,k ' )  has  to  be  synmet  lc:  H(k,k')  =  H(k,k'),  because 
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the  intersections  of  k  and  k'  are  identical  to  the 
intersections  of  k*  and  k. 

2.  In  order  to  minimize  the  number  of  comparisons,  we  use 
a  "gaod  guess"  of  which  patch  wi  .  be  visible  in  a 
given  interval.  This  guess  is  obtained  by  remembering 
the  "sample  points"  and  the  "sample  segments"  (see  also 
reference  [2])  of  the  previous  scan  line.  A  sample  seg¬ 
ment  is  a  segment  of  the  previous  scan  line  which  corres¬ 
ponds  to  a  visible  portion  of  a  conic.  A  sample  point 
is  the  extremity  of  a  sample  segment.  We  use  these 
points  of  the  previous  scan  line  in  oi der  to  get  sample 
points  belonging  to  the  same  patches  as  -  and  being  us 
close  as  possible  to  -  sample  points  of  the  previous 
scan  line. 

On  each  sample  interval,  we  have  to  compare  the  probably  visible 
patch  with  other  possibly  visible  patches  on  that  interval.  Thus  if 
there  are  n  patches  and  if  the  visible  segments  in  the  previous  scan 
plane  corresponded  to  v  different  patches  (v  <_  n) ,  then  the  number  of 
comparisons  will  be  less  than  (n  -l)v,  if  the  good  guesses  were  right, 
which  is  usually  the  case.  Moreover,  if  there  are  on  the  average  q 
possibly  visible  patches  behind  each  of  the  v  probably  visible  patches, 
then  the  total  number  of  comparisons  required  is  qv.  Thus  the  compute 
t  me  of  the  algorithm  is  proportional  to  the  "visible  complexity"  v  of 
the  picture. 

Of  course,  if  the  good  guess  is  wrong,  the  sample  segment  eventually 
has  to  be  divided  and  more  comparisons  are  needed: 

3.  If  A  is  in  front  of  B  and  3  is  in  front  of  C,  then  A 
is  in  front  of  C.  This  might  permit  us  to  avoid  having 
to  compute  the  intersection  of  A  and  C  in  some  cases 
where  the  intersections  of  A  and  B  and  of  B  and  C  were 
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already  computed  previously  on  the  scan  line.  The 
interest  of  this  "trick"  comes  from  the  fact  that  it 
is  much  more  expensive  to  compute  the  intersections 
of  two  conics  than  just  to  compute  their  depths  for 
some  value  of  x. 

4.  Another  improvement  which  was  suggested  to  the  author, 
but  nc  c  implemented  i  1  his  computer  program,  consists  of 
using  the  values  of  the  coordinates  of  the  intersections 
of  two  patches  in  the  previous  scan  plane,  as  initial 
approximations  of  their  values  m  the  current  scan  plane, 
and  to  refine  the  approximation  for  the  new  scan  plane  by 
applying  a  two-dimensional  Newton  method. 
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SHADING  THE  CURVED  SURFACES 

For  an  observer  looking  from  a  distance  d  in  the  z  direction,  and 
if  the  light  source  coincides  with  the  observer,  an  acceptable  and  pleasing 
shading  value  is  given  by: 

COS2CX 

(d-z)2 

At  the  point  of  the  surface  of  coordinates  (x,y,z) ,  a  is  the  angle 
between  the  normal  to  the  surface  and  the  direction  of  the  observer  (z  dir¬ 
ection)  }  we  have: 


cos  2  a  =  - - 

n2  +  n2  +  n2 
x  y  z 

where  (nx,  ny,  nj  are  the  components  of  the  normal,  given  by: 

"x  =  2alX  +  V  +  S5Z 

n  ■  2a.y  +  a.z  +  a,x 
y  2J  4  6 

n  =  2a, z  +  a,.x  +  a,y 
2  3  5  4 

Depending  on  the  type  of  quadric,  z  is  computed  either  by  solving 
a  quadratic  equation,  or  by  computing  the  ratio  of  the  value  of  a  second- 
order  polynomial  to  the  value  of  a  first-order  polynomial.  The  fact  that 
the  components  of  the  normal  and  the  value  of  z  have  to  be  computed  for 
hundreds  of  thousands  of  points  for  each  picture  makes  the  cost  of  shading 
generally  much  more  expensive  than  the  cost  of  visible  surface  detection. 
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THE  DATA  STRUCTURE 

An  interesting  by-product  of  this  investigation  is  a  system  „nich 
allows  one  to  create  and  modify  patches  ind  pictures  interactively. 

The  picture  designer  can  model  a  data  structure.  This  structure  can 
be  viewed  as  an  oriented  graph  without  loops,  in  which  appear  only  the  names 
of  objects  and  the  names  of  transformations,  and  not  the  actual  coefficients 
of  the  transformation  matrices  or  of  the  quadratic  equations.  The  same 
name  may  appear  in  several  places  in  the  structure,  and  these  occurrences 
do  not  have  to  be  changed  when  the  attributes  of  the  name  are  to  be  modified, 

The  nodes  of  the  graph  are  occupied  by  the  names  of  three  types  of 
objects : 

-  subpictures 

-  patches 

* 

an<3  -  quadratic  patches 

The  graph  edges  are  occupied  by  the  names  of  transformations  and, 
in  some  cases,  by  a  special  attribute  which  is  greater,  less  than,  or 
equ tl  to  zero.  This  special  attribute  appears  only  if  the  vertex  joins  the 
name  of  a  patch  to  the  name  of  a  quadratic  equation. 

A  language  has  been  designed  and  implemented  on  the  PDP-iO  to  allow 
one  to  define  transformations  and  quadric  patches,  and  to  create  or  delete 
objects  and  father-son  relationships.  Examples  of  sentences  of  this  lang¬ 
uage  are: 

total  AS  PROD  ROTX  0  CENTR  .5  .5  .5  TR2  .4  END 
which  defines  the  transformation  TOTAL  as  the  product  of  a  centered  rotation 
and  of  a  translation; 


MODE  NOINT 

QUA9R  CYL  AS  EQU  1  0  1  0  0  0  0  0  0  -1 
QliADR  CONE  AS  EQU  1  -1  1 
QUA  DR  PLANE  AS  EQU  0  0  0  0  0  (*  0  0  1 

TRANS  TCI  AS  PROD  SCALE  .1  TRX  .5  END 

TRANS  TC2  AS  PROD  SCY  2  ROTX  1S0  TR  .5 

Trans  tc?  as  prod  scale  .05  trx  .5  end 

TRANS  TC4  AS  PROD  SCY  2  ROTX  1F0  TR  .5 

TRANS  TBOTT  AS  ROTX  -90 
TRANS  TOPC1  AS  TL-OTT  CF.NTR  0  .5. 

TRANS  T0PC2  AS  TBOTT  CENTR  0  .6 

TRANS  TCPC3  AS  TLOTT  CENTR  0  .7 

TRANS  T0PC4  AS  TbOTT  CENTR  0  .8 

PATCH  Cl  IN  Cl  I NCL  CYL  BY  TCI 

I  NCL  PLANE  LY  TBOTT  PCS 
I  NCL  PLANE  BY  TCPC 1  NEC 
PATCH  CP.  IN  C2  I  NCL  CONE  LY  TC2 

I  NCL  PLANE  LY  TOPC1  POS 
I  NCL  PLANE  BY  TUPC2  NEC 
PATCH  C3  IN  C^  I  NCL  CYL  BY  TC3 

I  NCL  PLANE  BY  T0PC2  POS 
I  NCL  PLANE  LY  TCPC 3  NEC 
PATCH  C4  IN  C4  I  NCL  CONE  BY  TC4 

I  NCL  PLANE  BY  TOPC3  P  US 
I  NCL  PLANE  LY  T0PC4  NEC 
TRANS  OH  1  AS  PROP  ROTY  90  TRX  .5  END 
TRANS  OH 2  AS  PROD  CHI  R0T7.  223  CENTR  .3  .3  END- 
PATCH  A3  IN  A3  I  NCL  PLANE 

I  NCL  PLANE  LY  TBOTT  POS 
I  NCL  PLANE  LY  OH  1  POS 
I  NCL  PLANE  LY  OH2  POS 
I  NCL  CYL  BY  TCI  POS 
TRANS  TA2  AS  ROTY  -90  CENTR  .5 
SLIP  AT  IN  A?  I  NCL  A3  BY  TA2 

SILP  A!  IM  A I  I  NCL  A?  BY  TA2 

SUiP  A4  IN  -V  I  NCL  A1  LY  TA2 

SILP  TOTAL  IN  TOTAL  I  NCL  A1  I  NCL  A2 

I  NCL  A3  I  NCL  A  A 

I  NCL  Cl  I  NCL  C2 
I, NCL  C3  I  NCL  C4 
SCAN  50^ 

TRANS  TTT  AS  PROD  ROTX  30  ROTY  40  ROTZ  50  END  CENTR 
TRANS  TTT  AS  PROD  TTT  TH  -.1  .1  END 

MODE  SHADE 
SHADE  1  0  1200 

SO  P  ROTA  IN  ROTA  I  NCL  TOTAL  BY  TTT 
DIS.P  ROTA 
-TTY 


.7  END 
.8  END 


Figure  4 
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CONCLUSION 


Algorithms  have  been  presented  which  determine  the  visible  icrtions 
of  quadric  patches.  The  algorithm  for  patches  known  not  to  intersect  works 
much  faster  than  the  other  algorithm.  In  both  cases,  however,  the  computa¬ 
tion  of  the  precise  shading  value  at  each  point  of  the  picture  it  in  general 

the  most  expensive  process,  even  though  we  assumed  that  the  light  source  and 
the  observer  were  at  the  same  place. 

More  complicated  classes  of  surfaces  still  have  to  be  investigated, 
but  this  will  probably  require  more  computational  power.  It  would  also  be 
interesting  to  have  an  algorithm  which  (like  Watkins'  algorithm)  would  work 
very  fast  if  only  planar  polygons  are  involved,  but  which  could  also  handle 
quadric  surfaces  at  the  cost  of  a  smooth  degradation  of  performance. 
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